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1.  Introduction. 


A  suggestion  to  use  underground  cosmic-ray  measurements  for  the  detection  of 
in-situ  anomalies  of  rock  density  was  made  in  1976  by  Malmquist  et  al.  (Reference  2). 
The  method  is  based  on  the  dependency  of  cosmic-ray  intensity  on  the  average  rock 
density  above  the  detector.  Because  the  intensity  is  also  directional,  one  can  determine 
with  an  appropriate  detector  the  average  rock  density  in  a  specific  direction.  Malmquist 
et  al  proposed  to  use  the  method  for  geophysical  exploration,  particularly  for  the 
detection  of  ore  bodies.  Such  bodies  are  characterized  by  a  higher  density  than  the 
surrounding  rock  and,  therefore,  cause  anomalies  in  the  observed  cosmic-ray  intensities 
at  locations  that  are  deeper  than  the  ore.  The  method  has  also  been  used  for  the  search 
of  hidden  chambers  in  pyramids  (Alvarez  et  al.  1987,  Reference  3).  A  patent  for  a  bore¬ 
hole  apparatus  which  measures  directional  cosmic-ray  muon  intensities  has  been  issued 
in  1986  (Reference  4).  The  apparatus  could  provide  a  practical  instrument  for  the 
detection  of  small  anomalies,  such  as  cavities  (tunnels)  in  the  vicinity  of  a  bore  hole. 
Estimates  of  signal  strengths  which  can  be  expected  from  underground  tunnels  were 
presented  at  the  Third  Symposium  On  Tunnel  Detection  (Levy  et  al.  1988,  Reference  5). 
In  the  present  report  we  provide  more  elaborate  estimates  and  discuss  the  feasibility  of 
the  muon  radiation  method  for  the  detection  of  tunnels. 


2.  Normal  Radiation  Intensity. 

An  empirically  established  expression  for  normal  muon  radiation  intensity  in  a 
standard  Earth  is  (References  1  and  2) 


where 


and 


l(hp,0)  = 


I(hp.O)  =  l{hp,0)  cosn{hpk  , 

1.74-10® 


,-hp  8  10 


—7 


(AOO+hp  10-3)  (I1+/1/710"3)1  53 
n{hp)  =  1.85  +  4.65  10-6  (  hp\0~3  )*  67 


(2.1) 

(2.2) 

(2.3) 


In  these  equations, 

/  =  muon  counts  /  (  m2-s-  steradian  )  , 
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h  =  depth  below  the  ground  level  (Reference  2)  or  below  the  top  of  the 
atmosphere  (Reference  l),  m, 

p  =  average  density  of  matter  in  the  direction  of  measurement,  kg 'm3,  (the 
average  rock  density  at  moderate  depths  is  usually  assumed  to  be 
2700  kg  m3), 

8  =  zenith  angle,  radians. 

The  formulas  (2.1 )  and  (2.2)  are  taken  from  Reference  1.  They  have  been  established 
from  observations  of  radiation  intensity  up  to  depths  where  hp  =  7  106  kg/m2  The 
formula  (2.3)  for  the  exponent  of  the  cosine  is  derived  from  a  graph  of  the  function 
n(hp)  in  Figure  2  of  Reference  1.  The  formula  approximates  the  empirical  curve  for 
hp  <  2.7  10r  kg  m*  which  corresponds  to  about  1000  m  depth  For  depths  less  than 
200  m  the  exponent  is  practically  constant  and  equals  about  two 

The  complicated  dependence  of  the  normal  radiation  on  the  zenith  angle  6 
indicates  that  a  straight  ray  hypothesis  ("radiation  intensity  in  a  given  direction  is 
proportional  to  the  average  rock  density  in  that  direction")  is  only  approximately  true 
If  the  hypothesis  would  hold  exactly  then  at  moderate  depths  for  which  the  curvature  of 
the  surface  of  the  Earth  can  be  neglected  (say,  for  h  < 500  m  or  hp<1.35  106  kg  m2) 
one  would  have  a  radiation  formula  of  the  following  type: 

I(  hp,  6  )  =  I  ’  0 

cos  8 

Figure  1  displays  the  normal  radiation  intensity  as  computed  with  the  formulas  (2.1) 
and  (2.-1),  respectively.  Eq.  (2.1)  predicts  for  small  depths  a  smaller  radiation  intensity 
than  the  geometrical  formula  (2.-1)  whereas  for  large  depths  the  trends  are  reversed. 
The  reason  for  the  difference  at  large  depths  can  be  the  finite  radius  of  the  Earth  The 
reason  for  the  difference  at  small  depths  is  not  clear.  The  crossover  of  the  estimates 
occurs  at  about  60  metres 

For  the  present  calculations,  we  shall  assume  eq.  (2.1)  as  a  standard  and  use  it  for 
nonhomogeneous  rock  as  follows  (Reference  2).  First,  we  compute  the  average  rock 
density  in  the  direction  under  consideration,  and  then  use  that  density  and  the  vertical 
depth  h  in  the  formula.  For  simplicity,  we  shall  also  take  the  position  of  Reference  2 
and  assume  that  h  is  the  depth  below  ground  level  (i.e.  not  below  the  top  of  the 
atmosphere). 


3.  Effect  of  a  Tunnel  on  Observations  with  Zero  Aperture. 

We  assume  that  the  tunnel  is  an  underground  cavity  in  form  of  a  circular  cylinder 
with  its  axis  parallel  to  the  ground.  Let  the  point  of  origin  of  a  reference  coordinate 


system  be  located  at  the  radiation  detector  and  let  the  system  be  oriented  such  that  the 
'-axis  points  upwards  and  the  t/-axis  is  parallel  to  the  axis  of  the  cylinder.  In  this 
coordinate  system,  the  tunnel  (the  cylinder)  is  completely  described  by  its  radius  r0  and 
by  the  location  (.Y0  Z0)  of  its  axis  in  the  i, 2-plane  (Figure  2).  Let  the  7,2-plane  be  called 
the  reference  plane. 

Let  (t \6.6)  be  spherical  coordinates  related  to  the  Cartesian  coordinates  of  the 
reference  system  by 


7  =  v  ■  sin  8  ■  cos  6  , 
y  =  v  ■  sin  9  ■  sin  i>  , 
z  =  v  ■  cos  9 


The  spherical  coordinates  of  the  tunnel  axis  in  the  reference  plane  are 

v  =  R0  —\/ -Yfj  -f  Zl  . 

6  =0  . 

9  =  &0  =  arctan  (X0/Z0  )  . 


(3.2) 


The  support  half-angle  of  the  tunnel  cross  section  in  the  reference  plane  is 

J0  =  arcsin  f  r0  /  R0  |  .  (3.3) 

The  two  coordinates  <5  and  0  define  the  direction  of  a  ray  from  the  location  of  the 
detector  To  compute  the  effect  of  the  tunnel  on  radiation  received  from  that  direction 
one  needs  the  length  0(6.0)  of  the  tunel  chord  along  the  ray.  We  now  calculate  that 
length 


The  length  of  the  tunnel  chord  along  a  ray  (0,0)  in  the  reference  plane  is 


0(0.9)  = 


0 

2  (  r‘  -  /?0-  su r(9-60)  ) !/~ 


if  |0-©o!>^o 

if  1 0  —  0q  |  <  0O 


(3.-1) 


The  projection  of  a  ray  in  the  general  direction  (<2,0)  onto  the  reference  plane  has 
the  slope  cot  0/cos  6  ,  and  the  zenith  angle  of  the  projection  is 

90(6,9)  =  arctan  (  cos  <t>  tan  9  )  .  (3.5) 

The  projection  of  the  chord  D{6,9)  has  the  length  D(0, 90(6,9))  and  both  lengths  are 
related  by 
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(3.6) 


=  l/{u  ,  )  /  v  1  —  sin  o  sm  p  . 

Equations  (3.5)  and  (3.6)  allow  one  to  compute  the  lengths  of  tunnel  chords  along 
arbitrary  rays  from  the  point  of  origin. 

Let  ~rff(fr[Pr  be  the  depth  of  the  center  of  origin,  PTocjc  be  the  average  rock  density 
and  Ptunn; ;  be  the  density  of  matter  inside  the  tunnel.  Then  the  radiation  intensity  in 
the  direction  (6.0)  is  given  by  the  normal  radiation  intensity  formulas  (2.1)  through 
(2.3)  where  the  argument  hp  is  replaced  by 

det(  c‘  ''fP  rock  •  Ptunnel'  ^ zdelectorProck  (Ptunnel  Prock )  D(6,8)  COS  6  (30 

//  is  a  measure  for  the  mass  which  shields  the  detector  against  cosmic  radiation  in  the 
direction  (6.0).  Obviously,  if  the  rock  is  homogeneous,  that  is  P(unn(t=  P -ock-  then 
//  =  :ti(tcru>TPr0rk  The  absolute  effect  of  the  tunnel  on  the  radiation  intensity  is  the 
difference 

~  dttcct  orP  rock  '  Ptunnel-  *  -  6)  =  W  °)  ~  1  (uetectorProck •  *)  •  (38'> 

Its  dimension  is  counts  '(m2ssr):.  We  also  define  a  dimensionless  relative  effect  by 

E(:detcctorProck-  Ptnnnel'  0  ~  /  l(ZdetectorProck<  0)  ~  1  •  (3i^) 

The  last  argument  of  P  in  eq.  (3.9)  is  the  aperture  angle  of  the  detector  which  is 
assumed  to  be  zero  in  this  section. 

4.  Effect  of  a  Tunnel  on  Observations  with  Finite  Aperture. 

We  define  the  angular  aperture  of  the  detector  by  specifying  an  aperture  angle  o 
such  that  the  detected  radiation  is  received  from  within  a  cone  with  the  half-angle  a  (see 
Figure  3).  The  corresponding  solid  angle  of  the  aperture  is  27r(l  —  cos  a)  steradians. 
The  normal  radiation  which  is  received  with  the  aperture  a  equals  the  integral  of 
I ieieciorP rock  over  a^  ra)'s  (^at  are  ‘ns^e  the  aperture  cone.  The  effect  of  the 
tunnel  is  obtained  by  comparing  this  integral  with  the  integral  of  I(H,0)  over  the  same 
solid  angle. 

We  need  the  integration  limits  in  terms  of  the  spherical  coordinates.  They  are 
obtained  by  computing  the  intersection  of  the  aperture  cone  with  a  sphere  (with  its 
center  located  at  the  detector)  in  terms  of  the  coordinates  4>  and  6.  Let  tuc  radius  of  the 
sphere  be  R  and  the  axis  of  the  cone  be  given  by  the  coordinates  <t>c  and  6 f.  The  unit 
vector  in  the  direction  of  the  axis  is 
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sin  8C  •  cos  (t>c 
sin  8e  •  sin  d>c 
cos  8  e 


A  point  on  the  surface  of  the  sphere  is  defined  by  the  vector 


t?(W) 


R  ■  sin  9  ■  cos  <t> 
R  sin  8  sin  4> 
R  ■  cos  8 


(Al) 


(4.2) 


Points  of  intersection  of  the  sphere  with  the  aperture  cone  satisfy  the  equation 
(V  w)  =  R  cos  q  (see  Figure  3).  Pairs  (d,8)  which  define  rays  inside  the  aperture 
cone  sat isfy  the  inequality  fu  •!?)>/?■  cos  o .  or 

C(6.  8)  =  sin  6  sin  6C  ■  cos  (6  —  <5f)  +  cos  6  ■  cos  8 c  —  cos  a  >  0  .  (4.3) 

The  intensity  of  radiation  observed  with  the  aperture  a  in  the  presence  of  a  tunnel 
is 


Jl'detcctorProck'  Ptunnel'  *c'dc’a)  =  f  /(//,*)  '  sin  8  d*  d8  .  (4.4) 

C{o,0)  >  0 

J  has  the  dimension  [muon  counts  (m2s)[.  If  the  tunnel  is  not  present,  then  an 
observation  with  the  aperture  a  yields 

J^-deteclorProck'  Prock-<t,c-°e-a)=  f  zdete  ctorProck’9)  '  sin  6  d<t>  d6  '  5) 

C{<t>,e)  >  o 

The  absolute  effect  of  the  tunnel  is  for  the  aperture  a 
^ '(  *  dete  clorP  rock  •  P tunnel'  <*>'  9'  Q  )  — 

l  A 

—  ^  ‘’detector?  rock'  Ptunnel'  ^ '  9'  **)  9 ^detector?  rock  '  P  rock'  ^  ’  9'  ' 

The  dimension  of  F  is  jcounts/(m2s)[.  The  dimensionless  relative  effect  for  a  detector 
with  the  aperture  a  is 

F* '(  ‘‘detector? rock  '  Ptunnel’  <*>'  ° ) 

—  *^(  zdeteclor?rock  '  Ptunnel’  9’  **)  /  detector?  rock '  Prock’  ^  ^ 

Note  that  the  absolute  effect  F  depends  on  the  aperture  and  is  proportional  to  a2  if 
a«  1.  In  contrast,  the  absolute  effect  E,  eq.  (3.8),  is  for  a  single  ray  and  therefore  does 
not  depend  on  a. 
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5.  Examples. 


We  illustrate  the  effect  of  a  tunnel  on  the  observed  radiation  intensity  by 
calculating  the  radiation  for  a  number  of  typical  arrangements  of  tunnel  and  detector. 
To  be  specific,  we  assume  for  all  calculations  that  the  tunnel  diameter  is  2  m  and  that 
the  detector’s  aperture  a  is  5°.  The  normal  intensity  for  other  small  apertures  can  be 
estimated  from  the  presented  results  by  noting  that  the  intensity  is  proportional  to  a2. 
A  detector  with  a  large  aperture  produces  a  stronger  signal  than  a  detector  with  a  small 
aperture,  but  the  angular  resolution  of  the  device  with  the  smaller  aperture  is  better. 
This  means  that  the  smaller  the  aperture  the  better  is  the  accuracy  with  which  the 
tunnel  can  be  located.  Figu're  4  shows  that  an  aperture  of  about  5°  should  be  adequate 
for  depth  differences  up  to  20  m  between  tunnel  and  detector.  Much  larger  differences  of 
depths  would  require  smaller  apertures  in  order  to  achieve  a  reasonable  resolution. 
Figure  5  displays  the  normal  radiation  intensity  for  a  detector  with  an  aperture  of  5°. 
Note  that  the  intensity  in  Figure  5  has  a  dimension  different  from  that  in  Figure  1, 
because  Figure  5  shows  the  integral  (4.7)  of  the  normal  intensity  over  the  aperture  of 
the  receiver,  whereas  Figure  1  displays  the  normal  intensity  as  given  by  eq.  (2.1).  The 
general  trends  of  the  curves  in  both  figures  are  similar. 

If  the  aperture  is  finite,  then  the*  largest  effect  by  the  tunnel  likely  is  obtained  in  a 
direction  where  the  tunnel  is  closest  to  the  detector.  This  direction  is  in  the  plane 
orthogonal  to  the  axis  of  the  tunnel,  that  is  in  the  reference  x, e-plane  defined  by 
i>  =  0  mod(“).  We  test  this  property  of  the  tunnel  signature  by  calculating  the  intensity 
observed  in  directions  toward  points  of  the  tunnel  axis  at  a  distance  from  the  reference 
plane.  Figures  6  through  8  show  some  of  these  calculations.  They  were  made  for  a 
detector  depth  of  100  m,  a  tunnel  depth  of  80  m  and  two  tunnel  locations,  specified  by 
the  zenith  angles  6?0=  0°  and  45°.  In  the  case  of  zero  zenith  angle  the  tunnel  is  located 
directly  above  the  detector  at  which  position  it  has  the  largest  effect. 

Figure  6  shows  the  normal  intensities  observed  in  the  directions  of  the  two 
tunnels.  The  abscissa  is  the  distance  of  the  aim  point  from  the  reference  plane.  As  this 
distance  increases  (we  orient  the  detectors  aperture  field  towards  the  point  of  the  tunnel 
axis  indicated  by  the  abscissa)  the  corresponding  zenith  angle  increases,  too,  and 
consequently  the  normal  intensity  decreases.  The  effect  of  the  tunnel  cannot  be  seen  in 
the  plotting  scale  of  Figure  6  because  the  relative  effect  is  generally  less  than  1%  and 
therefore  mostly  within  the  line  thickness  of  the  curve.  The  relative  and  absolute  effects 
of  the  tunnel,  respectively,  are  shown  in  Figures  7  and  8.  The  relative  increase  of  the 
intensity  due  to  the  presence  of  the  tunnel  is  given  by  eqs.  (3.9)  and  (4.7),  and  the 
computed  results  are  shown  in  Figure  7.  The  dot-dash  curves  in  the  top  part  of 
Figure  7  show  the  relative  effect  of  the  tunnel  in  case  of  zero  aperture.  The  two  curves 
in  the  lower  part  of  the  figure  show  the  relative  effect  in  case  of  an  aperture  of  5°.  It  is 
interesting  to  note  that  for  zero  aperture  the  relative  effect  increases  slightly  as  the  axial 
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position  increases.  This  is  due  to  an  increase  of  the  tunnel  chord  for  oblique 
observations.  However,  this  effect  is  more  than  compensated  by  the  geometrical  effect  of 
the  finite  aperture  since  the  tunnel  occupies  a  smaller  part  of  the  field  of  view  as  the 
distance  to  the  tunnel  increases  The  absolute  tunnel  effect  is  shown  in  Figure  8  for  the 
same  two  tunnels.  The  figure  indicates  that  the  effect  of  the  tunnel  decreases  by  a  factor 
of  ten  as  the  distance  along  the  tunnel  axis  increases  to  about  40  m.  The  axial  distance 
of  40  m  approximately  corresponds  in  the  present  examples  to  a  zenith  angle  of  50°. 

In  summary.  Figures  6  through  8  confirm  that  the  largest  effect  of  the  tunnel  is 
indeed  in  the  reference  plane  for  all  zenith  angles  of  interest.  Also,  the  effect  is  reduced 
by  about  a  factor  of  ten  as  the  zenith  angle  increases  from  zero  to  about  45°.  Next  we 
present  some  calculations  of  the  tunnel  effect  in  the  reference  plane,  that  is  for 
conditions  optimal  for  tunnel  detection. 

Figure  9  shows  the  relative  effect  of  the  tunnel  in  the  reference  plane  for  three 

different  detector  depths  and  the  same  two  zenith  angles  of  the  tunnel  as  in  Figures  7 

and  8.  If  the  detector  is  placed  very  close  to  the  tunnel  then  the  outline  of  the  tunnel 
completely  covers  the  aperture’s  field  of  view.  The  corresponding  relative  effect  of  the 
tunnel  approximates  the  relative  effect  in  case  of  zero  aperture  which  is  independent  of 
the  tunnel  depth,  and  is  shown  in  Figure  9  by  dot-dash  lines.  In  case  of  the  5°  aperture, 

the  relative  effect  P  of  the  tunnel  is  reduced  approximately  by  a  factor  of  ten  when  the 

depth  difference  between  detector  and  tunnel  increases  from  zero  to  50  m. 

Figure  10  shows  the  absolute  effect  F  of  the  tunnel  for  the  same  cases  as  shown  in 
Figure  9  The  general  shape  of  the  curves  are  similar  to  those  of  Figure  9  except  that 
the  rate  of  decrease  of  the  absolute  affect  is  about  half  of  the  rate  of  the  relative  effect. 


6.  Discussion  of  Results. 

In  the  previous  section  we  presented  some  examples  of  the  effect  of  a  tunnel  on  the 
received  muon  radiation  intensity.  The  effect  was  found  to  be  of  the  order  of  one 
percent  and  the  accuracy  of  the  muon  counting  must  be  much  better  than  that  if  the 
tunnel  is  to  be  detected.  The  accuracy  of  the  muon  count  depends  on  the  total  count 
and,  since  the  total  count  depends  on  the  size  of  the  detector  and  duration  of  the 
observation,  one  can  use  calculations  of  the  effect  of  the  tunnel  to  find  estimates  of 
detector  sizes  and  observation  times  which  are  necessary  for  a  successful  tunnel 
detection  operation. 

The  strength  of  the  tunnel  signal  is  proportional  to  the  absolute  effect  F,  eq.  (4.6), 
illustrated  by  Figure  10.  In  order  to  detect  this  signal  it  is  necessary  that  the  standard 
error  of  the  observation  is  a  fraction  of  the  signal.  We  denote  the  fraction  by  e  and  call 
it  the  confidence  factor.  Let  the  standard  error  of  the  muon  count  be  aN,  and  let  the 
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increase  of  the  count  due  to  the  presence  of  the  tunnel  be  ANt unnt/-  Then  the  condition 
for  the  detection  of  the  signal  is 

oN  <  c  ANtunnel  .  (6  1) 


The  signal  strength  is 


*N,.nn«  -  F  I  -  ,  (6  2) 

where  t  [s]  is  the  duration  of  the  observation  and  a  [m2]  is  the  receiver  area  of  the 
detector.  The  standard  error  of  the  count  of  impacting  muons  equals  the  square  root  of 
the  total  count,  that  is 


oA  =  VJ  ■  t  ■  a  ,  (6.3) 

where  J  (counts  (m2s),  is  the  normal  radiation  intensity  for  the  given  aperture, 
eq.  (4.5).  Substituting  these  expressions  into  eq.  (6.1)  one  obtains  the  following 
condition  for  the  product  of  observation  time  and  detector  area: 

1  '  a  >  *7  "F  J  (M 

F- 


In  terms  of  the  relative  tunnel  effect  P  which  is  given  by  eq.  (4.7)  the  condition  for  the 
product  is 


t  a 


> 


J  P 2 


F  P 


(6.5) 


A  reasonable  value  of  the  factor  \  A2  is  four,  corresponding  to  t  =0.5  or  to  a 
signal  that  is  twice  ar  large  as  standard  error  of  muon  counting.  A  more  reasonable 
value  of  the  factor  is  ten,  that  is  €  =  0.316  which  roughly  means  that  the  signal  is  about 
three  times  as  large  as  the  standard  error  of  counting.  In  our  examples  we  shall  use  the 
value  ten  for  l/c2.  Therefore,  the  estimated  conditions  for  the  product  t  a  might  be 
relaxed  possibly  by  a  factor  of  two,  although  this  is  not  recommended.  (We  have  only 
considered  the  counting  accuracy.  Additional  measurement  noise  can  be  caused  by 
temporal  and  spacial  variations  of  the  muon  radiation  as  it  arrives  at  the  Earth,  and  by 
changing  conditions  in  the  atmosphere  during  the  recording). 


The  product  t  a  is  a  measure  of  the  effort  that  is  needed  to  produce  the  required 
result.  Let  it  be  called  observation  effort.  Its  dimension  is  s  m2,  which  we  replace  by  the 
more  practical  dimension  days  m2.  Next,  we  compute  the  required  observation  effort  for 
the  cases  shown  in  Figures  9  and  10. 


Figure  11  shows  the  required  observation  effort  for  detector  depths  of  30,  50,  100 
and  200  m.  To  determine  whether  the  required  effort  for  these  cases  has  a  manageable 
order  of  magnitude,  one  has  to  estimate  the  detector’s  receiving  area  a.  For  instance,  if 
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the  area  a  is  about  one  square  metre  then  a  tunnel  detection  with  the  receiver  at  100  m 
depth  would  require  observations  which  last  months  or  even  years.  To  simplify  the 
discussion  we  assume  that  the  receiving  area  equals  the  silhouette  of  the  detector,  and 
estimate  the  measurements  of  such  a  device.  (Actually,  the  effective  receiving  area  could 
be  a  fraction  of  the  silhouette,  or  a  combination  of  several  silhouettes,  depending  on  the 
construction  of  the  device  and  on  the  data  recording  logic). 

The  properties  which  distinguish  a  tunnel  from  geological  density  variations  are 
well  defined  boundaries  and  the  magnitude  of  the  density  anomaly.  If  the  extension  of 
the  device  is  much  larger  than  the  tunnel  diameter  then  the  signature  of  the  tunnel  is 
diminished  and  smeared  out,  and  it  might  not  be  possible  to  distinguish  between 
geological  density  variations,  topographic  effects  and  the  tunnel  in  the  tomographically 
reconstructed  density  field  Therefore,  we  require  that  the  linear  extension  of  the 
detectors  silhouette  is  not  more  than  2  m  (the  tunnel  diameter)  in  a  direction  orthogonal 
to  the  tunnel  axis.  Assuming  that  the  zenith  angle  of  the  observations  is  less  than  45° 
one  obtains  that  the  vertical  dimension  of  the  detector  should  be  less  than  2/sin  45°  or 
2.82  m  The  radial  dimension  of  the  detector  is  bounded  by  the  diameter  of  the  bore 
hole  in  which  the  detector  is  placed 

The  detector  proposed  in  Reference  4  has  the  form -of  a  cylinder  aligned  with  the 
axis  of  a  bore  hole  The  device  registers  all  those  muons  and  the  directions  of  their 
paths  which  penetrate  the  side  walls  of  the  cylinder.  Thus  one  can  assume  that  for 
vertical  bore  holes  the  apparatus  will  have  a  directional  sensitivity  proportional  to  sin  9C 
where  9 (  is  the  zenith  angle  of  the  direction  of  the  incoming  muons.  Let  the  length  of 
the  detector  be  L  and  its  diameter  be  D.  Then  the  effective  recording  area  is 
approximately  given  by 

a  =  D  ■  L  •  sin  6C  .  (6.6) 

Following  Reference  5  we  assume  that  the  bore-hole  diameter  is  0.30  m  and  the 

corresponding  effective  diameter  of  the  detector  is  D  —  0.26  m.  From  the  resolution 

requirement  discussed  above  we  have  L  =  2/sin  9C  m3X,  where  the  maximum  zenith 
angle  for  the  observations  might  be  about  45°.  Thus,  the  effective  recording  area  is 
given  by 

a  =0.26  ,  ■  sin  9  <  0.52  [  m2  1  .  (6.7) 

Sln  6c  max 

One  sees  from  Figure  11  that  with  this  recording  area  the  required  recording  time  is 
tens  or  even  hundreds  of  days  for  detector  depths  of  50  m  or  more.  Hence,  the  proposed 
tunnel  detection  method  and  device  cannot  be  reasonably  used  for  depths  over  50  m. 
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Figure  12  shows  plots  of  required  observation  effort  for  some  detector  depths  less 
than  50  m.  The  plots  indicate  that  a  practical  bound  for  detector  depths  is  between  30 
and  40  m.  The  required  recording  times  at  those  depths  are  of  the  order  of  days  or  tens 
of  days. 

To  determine  the  bound  more  accurately  we  assume  a  specific  receiver  area  and 
compute  the  corresponding  required  recording  times  in  days  for  detector  depths  of  30  m 
and  40  m,  respectively.  To  compute  the  receiver  area  we  use  eq.  (6.7)  with  the 
maximum  zenith  angle  of  45°.  (Thus,  the  results  are  for  a  bore-hole  diameter  of  30  cm 
and  a  detector  length  L  =  2.82  m).  Figure  13  shows  the  results  for  a  30  m  detector 
depth  and  tunnel  depths  between  5  m  and  25  m.  Because  we  must  assume  that  the 
tunnel  depth  is  not  known,  the  recording  times  have  to  be  such  that  the  worst  case  is 
covered.  In  the  present  example,  this  is  the  case  where  the  tunnel  is  at  a  depth  of  5  m, 
that  is  the  top  curve  in  Figure  13.  (Extrapolation  to  a  shallower  tunnel  at,  e.g.  1.5  m 
should  be  obvious  from  the  figure).  The  curve  shows  that  one  needs  about  ten  days  of 
observation  time  if  the  observation  zenith  angles  are  restricted  to  the  interval  between 
5°  and  45°.  This  result  indicates  that  30  m  is  a  practical  upper  bound  of  the  detector 
depth  for  tunnel  detection  operations.  The  corresponding  maximal  depth  of  the  tunnel 
depends  on  the  distance  between  the  bore  holes  from  which  the  measurements  are  used 
for  the  tomographic  density  reconstruction.  In  any  case,  the  tunnel  depth  must  be  less 
than  the  depth  of  the  detector. 

To  see  how  much  the  requirements  change  when  the  detector  is  placed  at  greater 
depths  we  plot  in  Figure  14  the  required  observation  times  for  a  detector  depth  of  40  m 
and  tunnel  depths  between  5  m  and  35  m.  The  worst  case  requirement  (read  from  the 
5  m  curve)  is  now  about  40  days  of  observation  time,  which  seems  to  be  too  long  for 
practical  operations. 

The  presented  estimates  of  required  observation  times  depend  on  a  number  of 
assumptions  about  the  mode  of  operation  and  properties  of  the  detector.  We  now 
discuss  consequences  of  differing  assumptions  about  the  tunnel  detection  method.  First, 
a  reduction  of  observation  times  by  about  a  factor  of  two  could  be  achieved  by 
restricting  the  zenith  angle  interval  of  the  observations  because  the  signal  of  the  tunnel 
is  stronger  for  smaller  zenith  angles  (see  Figures  9  -  12).  A  consequence  of  such  a 
restriction  is  a  reduction  of  the  explored  area  which  means  that  more  bore  holes  at 
closer  distances  would  be  needed  to  cover  the  same  area.  Second,  the  observation  times 
could  be  reduced  if  the  observations  were  made  with  larger  apertures  and  longer 
detectors.  The  consequence  of  an  increase  of  the  aperture  and  detector  length  is  a 
reduction  of  spacial  resolution,  that  is  an  impairment  of  the  detectability  of  the  tunnel 
by  density  reconstruction.  Third,  the  confidence  factor  <  could  be  increased  from  0.316 
to,  say  0.5.  As  discussed  above,  this  would  decrease  the  required  observation  times  by 
about  a  factor  of  two.  Because  this  reduction  of  observation  times  reduces  the  signal  to 
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noise  ratio,  it  also  reduces  the  probability  of  a  detection  of  the  tunnel.  Finally,  a  small 
reduction  of  the  observation  time  could  be  achieved  by  designing  the  detector  such  that 
its  sensitivity  is  not  proportional  to  the  sine  of  the  zenith  angle  thereby  enhancing  the 
sensitivity  at  small  zenith  angles.  It  seems,  however  that  any  combination  of  these 
modifications  could  possibly  reduce  the  required  observation  times  by  not  more  than  an 
order  of  magnitude  which  would  increase  the  practical  bound  for  the  detector  depth 
from  30  in  to  only  about  40  m. 


The  conclusions  in  this  report  differ  significantly  from  the  conclusions  by  Levy  et 
al  in  Reference  5.  Possible  reasons  for  the  differences  are  different  formulas  for  the 
normal  radiation,  approximate  calculations,  different  detector  geometries,  and  different 
definitions  of  the  effect  of  the  tunnel.  We  now  discuss  in  turn  these  possible  sources  of 
disagreement. 


The  normal  radiation  formula  used  by  Levy  et  al.  is 


I(hp,9)  =  I 


cos  9 


0 


1 

cos  9 


(6.8) 


This  formula  differs  from  the  formula  (2.1)  used  in  this  report  and  from  the  geometrical 
formula  (2.4),  and  it  is  quoted  by  Levy  et  al.  without  derivation  or  reference.  From 
Figure  1  where  the  geometrical  formula  is  compared  with  eq.  (2.1)  (for  6  =  45°)  one  can 
conclude  that  eq  (6.8)  predicts  for  shallow  detector  depths  up  to  three  times  larger 
radiation  intensities  than  eq.  (2.1).  However,  the  difference  between  the  formulas  is 
insignificant  for  a  detector  at  100  to  200  m  depth  which  is  the  principal  range  of  interest 
in  Reference  5. 


To  simplify  our  calculations  we  assumed  that  the  detector  length  is  small 
compared  with  the  distance  to  the  tunnel.  Levy  et  al.  do  not  make  that  assumption  but 
do  make  a  number  of  other  simplifying  assumptions.  The  precise  effect  of  these 
assumptions  is  difficult  to  assess  and  a  direct  comparison  of  the  numerical  results  is  not 
possible  because  of  differing  assumptions  about  the  detector  size  and  definitions  of  the 
tunnel  effect.  An  order  of  magnitude  comparison  is,  however,  possible  and  it  is  shown  in 
the  Appendix  that  results  in  both  papers  agree.  Hence  one  might  conclude  that  the 
various  simplifying  assumptions  in  both  papers  have  not  distorted  the  results. 

We  assumed  in  this  paper  that  the  vertical  length  of  the  detector  is  2.82  m  and 
that  the  radiation  is  recorded  with  a  5°  aperture.  We  then  calculated  the  radiation 
within  the  corresponding  solid  angle  with  and  without  the  tunnel,  and  defined  the  tunnel 
effect  as  the  difference  between  both  computations.  (The  calculations  were  done  in  the 
reference  plane,  that  is  for  the  azimuth  where  the  tunnel  has  the  largest  effect).  Levy  et 
al.  assume  that  the  detector  extends  from  the  ground  surface  to  200  m  depth.  Instead  of 
defining  an  aperture  cone  they  estimate  the  solid  angle  which  is  supported  by  the  tunnel 
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(shown  as  a  horizontally  hatched  area  in  Figure  4),  calculate  the  radiation  within  this 
angle  with  and  without  the  tunnel,  and  define  the  tunnel  effect  as  the  difference  between 
these  two  calculations.  Because  of  this  definition,  their  relative  as  well  as  absolute 
tunnel  effects  are  much  larger  than  ours  and,  consequently,  their  estimates  of  required 
observation  times  are  much  shorter. 

It  is  not  obvious  how  the  calculations  by  Levy  et  al.  are  related  to  measurements 
which  are  needed  for  a  tomographic  density  reconstruction.  In  general,  one  needs  for 
such  a  reconstruction  the  angular  dependence  of  the  signal,  and  observations  from  a 
number  of  receivers  at  different  locations.  These  requirements  dictated  our  choice  of  5° 
for  the  aperture  in  order  to  register  the  angular  dependence  of  the  signal,  and  2.82  m  for 
the  detector  length  in  order  to  have  a  well  defined  detector  location.  (See  the  beginning 
of  Section  5  and  the  discussion  of  eq  (6.7)  ).  Measurements  with  finite  apertures  for 
tomographic  reconstruction  are  also  suggested  by  Levy  and  Mockett  in  Reference  4.  In 
order  to  realize  the  tunnel  effect  estimated  by  Levy  et  al.  in  Reference  5  one  must  either 
know  the  tunnel  location  or  register  the  radiation  in  a  very  wide  angle  of  aperture. 
Wide  aperture  observations  not  only  smear  out  the  angular  dependency  but  also 
increase  the  total  count  and  the  corresponding  required  observation  times  (see  eq.  (6.4)  ) 
over  the  estimates  provided  by  Levy  et  al.  Also,  the  nominal  depth,  that  is  the  location 
of  the  observer  in  the  reference  plane  is  not  very  well  defined  if  the  detector  extends 
from  the  ground  surface  to  a  depth  of  200  m.  This  uncertainty  possibly  can  reduce  the 
accuracy  of  a  reconstructed  density  distribution. 

7.  Conclusions. 

Subterranean  cosmic-ray  intensity  depends  mainly  on  the  average  density  of  the 
material  which  has  been  penetrated  by  muons.  Therefore,  measurements  of  muon 
radiation  intensity  can  be  used  to  locate  density  anomalies  in  the  overlaying  strata. 
Such  anomalies  can  represent  either  geological  structures  or  man  made  artifacts,  like 
tunnels.  This  method  of  locating  density  anomalies  has  been  proposed  in  Reference  5 
for  the  detection  of  tunnels,  whereby  the  detector  apparatus  would  be  located  in  a  bore 
hole.  The  measurements  from  a  multitude  of  detectors  would  be  combined  and  the 
density  distribution  above  the  detectors  determined  by  tomographic  reconstruction. 
The  feasibility  of  this  method  for  tunnel  detection  depends  on  the  recording  times  which 
are  required  for  the  measurements.  In  this  report,  the  effects  of  cylindrical  tunnels  with 
diameters  of  2  m  on  the  radiation  intensity  are  calculated  and  the  required 
measurement  efforts  estimated.  The  conclusion  from  these  estimates  is  that  the  method 
might  be  feasible  for  tunnel  detection  up  to  detector  depths  of  not  more  than  40  m.  At 
these  depths,  the  necessary  recording  times  are  of  the  order  of  days,  approaching  the 
order  of  a  month  in  limit  cases.  If  the  depth  of  the  detector  is  increased  then  the 
corresponding  required  recording  times  increase  exponentially.  The  depths  of  detectable 
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tunnels  are  less  than  the  depths  of  the  detectors. 

The  average  rock  density  was  assumed  to  be  2700  kg/m3  for  the  computation  of 
tunnel  effect  estimates.  Deviations  from  this  average  by  some  10%  are  possible  under 
normal  circumstances,  and  they  would  change  the  depth  estimates  proportionally. 
Because  of  the  approximate  nature  of  the  estimates,  changes  of  this  order  of  magnitude 
can  be  considered  as  negligible. 

The  estimates  are  based  on  simple  assumptions  about  the  properties  of  the 
radiation  detector  and  requirements  on  data  for  a  tomographic  density  reconstruction. 
It  is  conceivable  that  special  data  analysis  methods  can  be  developed  which  use  prior 
knowdedge  about  the  tunnel  (for  instance,  its  direction  and  size)  and  are  less  demanding 
on  data  accuracy  than  a  general  tomographic  density  reconstruction.  Therefore,  a 
decision  to  support  a  construction  of  the  proposed  detector  for  application  in  tunnel 
detection  should  be  based  on  results  from  a  tested  data  interpretation  program  on 
simulated  data.  Such  results  would  not  only  establish  more  accurate  bounds  of 
applicability  of  the  method  but  also  provide  design  parameters  for  the  detector. 
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Appendix.  Comparison  of  Numerical  Results. 

We  compare  numerical  results  by  Levy  et  al.  in  Reference  5  with  this  report  to 
check  the  calculations  and  the  adequacy  of  simplifying  assumptions. 

In  our  calculations,  the  detector  diameter  is  0.26  m  and  the  detector  length  is 
2.82  m.  The  dimensions  of  the  detector  are  assumed  to  be  short  compared  to  the 
distance  to  the  tunnel.  The  effect  of  the  tunnel  we  define  as  the  difference  between 
observations  with  and  without  the  tunnel  within  a  solid  angle  corresponding  to  a  cone 
with  5°  half  angle.  A  summary  of  the  required  observation  effort  is  shown  in  Figure  11. 


Levy  et  al.  assume  in  their  numerical  example  a  detector  with  a  diameter  of 
0.26  m  and  a  length  that  extends  from  the  ground  surface  to  200  m  depth.  They  define 
the  tunnel  effect  as  the  difference  between  observations  with  and  without  tunnel  within 
the  solid  angle  that  is  supported  by  the  silhouette  of  the  tunnel.  A  quantitative 
summary  of  their  results  is  shown  in  Figure  2  of  Reference  5  which  displays  those 
locations  of  tunnels  where  the  effects  are  sufficiently  large  to  be  detected  in  seven  days. 

To  compare  the  results  one  should  choose  comparable  geometries.  Therefore  we 
select  a  tunnel  that  is  close  to  the  detector  so  that  the  5°  half  angle  of  our  aperture  cone 
approximately  corresponds  to  the  support  half  angle  of  the  tunnel  cross  section.  The 
tunnel  depth  we  choose  close  to  200  m  so  that  the  portion  of  the  200  m  long  detector 
which  registers  its  effect,  i.e.  the  portion  that  is  deeper  than  the  tunnel,  is  short. 


Specifically,  we  assume  Z0=  190  m  and  A'0=5m,  and  read  from  Figure  2  of 
Reference  5  that  the  tunnel  effect  equals  about  1.8  for  an  observation  length  of  seven 
days.  Because  only  those  detector  parts  which  are  deeper  than  the  tunnel  register  a 
tunnel  effect  we  estimate  the  effective  depth  of  the  detector  to  about  200  m,  and  its 
effective  length  to  about  5  m.  Then  the  corresponding  recording  area  is  about 
5  0. 26  =  1. 3  m2.  The  support  half-angle  of  the  tunnel  cross  section  is  about  arctan(l/10). 
Levy  et  al.  estimate  that  most  of  the  tunnel  effect  is  observed  within  a  solid  angle  that  is 
supported  by  the  tunnel  silhouette  with  a  length  twice  the  nearest  distance  to  the  tunnel 
(p.  288).  The  size  of  this  solid  angle  which  they  use  for  their  radiation  estimates  is 


Al  « 


2  ■  arctan  — 
10 


|  2  ■  arctan  1  j 


0.31  sr 


(Al) 


For  2je(ector  =  200,  Z0  =  190  and  X0  =  5  we  obtain  60  =  26°.  Figure  11  provides 
for  this  input  a  required  observation  effort  of  about  300  [days  m2j  if  the  tunnel  effect 
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should  be  larger  then  3.16  aN.  The  solid  angle  of  the  observational  field  is  for  Figure  11 
assumed  to  be 


Ac  =  2  jt  (1  —  cos  5°  )  =  0.024  sr  . 


(A.2) 


Scaling  our  result  from  Figure  11  to  the  conditions  in  the  Levy  et  al.  paper  we  obtain 
for  the  required  number  of  days  the  value 


!\r  «  300  jdaysm2]  • 


1-8 

3.16  cJm 


0.024  jsrl 
0.31  Isrl 


1.3  [m2l 


6  days 


instead  of  N  —  7.  The  estimate  (A. 3)  can  be  increased  or  decreased  by  a  large  factor 
because  the  estimates  of  the  effective  detector  length  and  aperture  solid  angles  arguably 
can  be  different  and  the  reading  of  1.8  a ^  from  Figure  2  of  Reference  5  is  quite 
inaccurate.  Therefore,  it  is  easy  to  obtain  an  exact  agreement  between  both  papers.  We 
conclude  from  the  comparison  that  the  results  in  both  papers  agree  in  order  of 
magnitude. 
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Figure  1.  Normal  radiation  intensity. 
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Figure  2.  Tunnel  coordinates. 
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Figure  3  Aperture  angle  o c 
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Figure  5.  Normal  radiation  intensity 

for  a  finite  aperture. 
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Axial  position  in  tunnel,  m 


Rock  density  2700  kg/m3 
Detector  depth  ■■  1 00  m 
Detector  aperture  »  5° 
Tunnel  depth  «  80  m 
Tunnel  radius  —  1 .0  m 


Figure  6.  Normal  radiation  intensity 

in  tunnel  direction. 
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Figure  7.  Relative  tunnel  effect 

versus  axial  position. 
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Figure  8.  Absolute  tunnel  effect 

versus  axial  position. 
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Figure  9.  Relative  tunnel  effect 

in  th.e  reference  plane. 
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Figure  10.  Absolute  tunnel  effect 

in  the  reference  plane. 
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-  ■—  zenith  angle  zero 

Confidence  factor  cr/s ignal  —  0.316 

Figure  11.  Required  observation  effort 

in  the  reference  plane. 
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Figure  12.  Required  observation  effort 

in  the  reference  plane 
for  shallow  tunnels. 
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Rock  density  ■■  2700  kg/m3 
Detector  aperture  —  5° 

Detector  depth  ■■  30  m 

Tunnel  radius  ■»  t.O  m 
Confidence  factor  a/stgnal  ■■  0.316 
Bore  hole  diameter  ■■  30  cm 

Figure  13.  Required  observation  time 

in  the  reference  plane 
for  a  detector  at  30  xn 
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Figure  14.  Required,  observation  time 

in  the  reference  plane 
for  a  detector  at  40  zn 
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